How delicate are the f(R) gravity models with disappearing 

cosmological constant? 
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We consider stability of spherically symmetric solutions in f(R) gravity model proposed by 
Starobinsky. We find that the model suffers from a severe fine tuning problem when applied to 
compact objects like neutron stars. The problem can be remedied by introducing a cut off on the 
mass of the scalar degree of freedom present in the model. A new mass scale associated with neutron 
stars density is then required for the stabilities of f(R) gravity solutions inside relativistic stars. 
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I. INTRODUCTION 



The cause of cosmic repulsion responsible for late time 
acceleration is one of the mysteries of modern cosmol- 
ogy. The simplest possibility to account for this effect is 
related to the assumption of dark energy Q. It is. how- 
ever, quite possible that late time acceleration of uni- 
verse is the result of large scale modification of gravity. 
Amongst all the schemes of modification of gravity in the 
infrared regime, the f(R) theories of gravity Qare most 
elegant and promising. These theories apart from the 
spin-2 object necessarily contain a scalar degree of free- 
dom dubbed scalaron. Stability of the theory requires 
that the scalaron is not tachyon and graviton is not a 
ghost which can be ensured by demanding the positiv- 
ity of the first and the second derivatives of f(R) with 
respect to the Ricci scalar R. 

The local gravity constraints impose most stringent re- 
strictions on any scheme of large scale modification of 
gravity in particular the f(R) theories of gravity. Most 
of the f(R) models are either not cosmologically viable or 
simply reduce to ACDM0. An interesting class of mod- 
els proposed by Hu-Sawicki and Starobinsky , referred to 
as HSS hereafter, can reconcile with local gravity con- 
straints and has the potential capability of being distin- 
guished from cosmological constant 0,1a @ (See Ref.0] on 
the related theme). 

The de-Sitter minimum in these models is very near 
to the curvature singularity which the scalaron can 
easily hit while evolving towards the minimum of its 
potential!!, |, El EH EL El EL E|- 0ne can tr y 
to modify HSS models by creating a large potential 
barrier between the de-Sitter minimum and curvature 
singularity (l6|. However, the latter results arc in clear 
violation of local gravity constraints [l7j- Thus the pres- 
ence of finite time singularity is generic in these models 
and should be handled carefully. The safe passage of 
scalaron to the minimum of its potential requires fine 
tuning of its initial conditions. The situation gets worse 
in the high curvature regime. 

In this paper we examine the fine tuning problem asso- 
ciated with f(R) theories with disappearing cosmological 
constant. We demonstrate that the level of fine tuning 
the HSS models require in case of the compact objects 



like neutron stars poses serious challenge to these mod- 
els. The problem can be alleviated by introducing the 
quadratic curvature terms in the HSS scenario. 



II. THE SCALAR DEGREE OF FREEDOM IN 
f(R) THEORIES OF GRAVITY 

We consider the modification of the Einstein-Hilbert 
action in the presence of the matter Lagrangian C m 



S 



(1) 



where f(R) is a function of the Ricci scalar R. Varia- 
tion with respect to metric leads to the following field 
equations 



□/,« - If) 9„ v = T^, (2) 



where f tR = df/dR and T^„ = -25C m /5g^ + g^Cm- 
The f(R) theories of gravity necessarily contain an addi- 
tional scalar degree of freedom which becomes clear after 
taking the trace of 10 



dV 



(3) 



where $ = f M ; dV/dijj = (2/(i?(V0) - ^R^)) /3. 

In what follows, it would be convenient to us to write 
f(R) in the form of the correction term to Einstein- 
Hilbert action, A, 



f(R)=R + A, 



1 + A. 



(4) 



where A ^ denotes the derivative of the correction term 
with respect to the Ricci scalar R. 

In case of Starobinsky model 0] , we have 



A = -XR 



( § 



1 - (1 + R 2 /R 

2ra+l 



(5) 

for R > R c . (6) 



Taking R c = p A /A where p A is the cosmological density, 
(10 -29 g/cm 3 ) and R ~ p c is the curvature inside the 
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neutron star such that (p c ~ 10 14 g/cm 3 ). The numerical 
value of ipo corresponding to the de-Sitter minimum is 
given by 



1-O(l0- 122 ) n = 0.9, 
1-O(l0~ 217 ) n = 2. 



(7) 



The Local Gravity Constraints are satisfied for 71 > 0.9 
fl8l ] and the evolution of density perturbations during 
the matter-dominated epoch requires n > 2 for G e // = 
G/f'(R) to be consistent with observations [f|. 

The scalar degree of freedom plays an important role 
in f(R) theories of gravity namely, its dynamics controls 
the space-time curvature. In generic cases, the de Sitter 
minimum at ip® is very close to ip = 1 corresponding 
to curvature singularity. The finite barrier between the 
singularity and de Sitter minimum means the curvature 
singularity is energetically accessible. In order to avoid it, 
the evolution of the field needs to be fine tuned. We shall 
see later that in the case of the neutron star with constant 
density p c , the extreme fine tuning of initial conditions 
becomes necessary for the existence of GR-like solution 
(R ~ p c ) along the radius r of the star to match the 
correct boundary conditions at its surface, r = r*. 



III. THE GROWING MODE OF THE 
PERTURBATION AND FINE TUNING OF 
INITIAL CONDITIONS 

The problem of fine tuning is inherent in /(i?) theo- 
ries if they are to be consistent with local gravity con- 
straints and can be appreciated by using the analytical 
arguments. An approximation scheme for the solution of 
ip(r) can be set up using the iterative computation for 
R(r) as followsjJlH, 



R{r) = Ro(r) + MZi(r) + 8R 2 {r) + 



(8) 



where Ro(r) — p(r) — 3P(r) and the first order iteration 
gives rise to the following expression, 

SR X = [-3V 2 A + A,rR — 2A] | r=Rq . (9) 

In case of Starobinsky model, we have 



|V 2 A,*|,|A, fl i?|«|A| 
«?i(r) w -2 A fu 2XR C 



for R > R c , (10) 
const. (11) 



In the first order iteration, ipi can be expressed through 
the GR-like solution tpo as 



ipi 



2n+l 



2nX 



^ +4nA 2 (2n+l) 



[I + SRx/Rq] 

2n+2 



-2n-l 



(12) 
(13) 

(14) 



where we have used 5R\ = 2XR C . 

Let us note that the scalaron mass in high curvature 
regime R « Rq ~ p ^> R c is given by 



fv, 

dip' 



R, 



\4>=4>o 



6n(2n+ 1)A \R, 



2n+2 



, (15) 



which allows us to obtain the first order iteration solu- 
tion, 



2 R 
6 m 



(16) 



The GR-like solution, ipi(r) under consideration, can de- 
viate from ipo(r) only near the stellar radius otherwise 
many known observational constraints of neutron stars 
will not be satisfied. Indeed, the first order iteration 
ipi solution is approximately the Schwarzschild de Sit- 
ter solution because it is corresponding to the curvature 
R — 2XR C which differs from R of the GR solution by a 
constant. In large scalaron mass limit, ip reduces to GR 
solution as expected. 

The configuration of the perturbative solutions of the 
f(R) gravity and general relativity are very different in 
the limit of the large scalaron mass. As demonstrated 
Refs.fl9| fc[20l|. the sign and the size of play a crucial 
role for the stability of solutions in time. The Dolgov- 
Kawasaki instability can be avoided by choosing m 2 , > 
which makes the perturbation Sip(t) oscillating in time. 
However, a large positive m 2 causes the instability of 
static spherically symmetric solutions [2l| for a class of 
f(R) gravity models which are carefully built to evade 
the local gravity constraints Q 1 ■ 

To demonstrate the catastrophic instability from a 
large positive value m 2 , let us consider a small pertur- 
bation Sipi(r) around ^p\(r). Assuming the static spher- 
ically symmetric metric, the trace equation ([3]) tells us 
that 



5ip" + -5ip[ 
r 



dV 
# 
d 2 V 



dV 

^=ip 1 +5ip 1 dip 



Sipi 



(17) 
(18) 



where primes denote derivatives with respect to r. 

In case m 2 ^ > 0, the growth of the perturbation Sipi (r) 
along the radius can be obtained in form 



5Mr) = SMO) \Ci + C 2 

r r 



(19) 



where we have used the notation mo = m^, 1 m m^ . 
It follows from Eq. (Tl9|) that the exponentially growing 



1 This type of instability that we shall focus on in the subsequent 
discussion refers to the radial evolution. 
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mode of 6ip% is unavoidable. Thus the instability of solu- 
tions always persists for any metric f(R) gravity models 
with the large m^. This fact may be exhibited by cast- 
ing the equation of perturbation 6ipi(t, r) in the following 
form, 



(df - V 2 ) Sip = -m%5il>. 



(20) 



It should be noticed that the difference of the sign of <9 2 
and V 2 comes from the signature of the metric itself. The 
avoidance of Dolgov-Kawasaki time instability by choos- 
ing m 2 , > invokes the instability (Sipi(r) oc e ±m ° r ) of 
the static solution. The orthogonality of the stability 
conditions of the time and space arising from the sig- 
nature of the metric was proposed in Ref. [2l[ and the 
evidence of this instability of the static solutions was ob- 
served numerically as the problem of the existence of rel- 
ativistic stars in f(R) gravity theories [H|,[l2|- 

It is however difficult but necessary to maintain the 
small deviation from the GR-like solution ipo(r) for the 
whole range of the stellar radius when the growth of per- 
turbation is exponential. The initial value 6ipi(0) must 
be extremely fine tuned if we want to stay near the GR- 
like solution. The seriousness of the fine tuning is re- 
lated to the size of the number mor* which is typically 
huge. The length scale corresponding to mo, the Comp- 
ton wavelength A c = I/thq, is very small in the nuclear 
matter density regime due to chameleon effect. Without 
the cut off on mo, A c can shrink below Planck length( 
r p = y/ KG I c 3 ~ 10 -33 cm). Meanwhile, the stellar ra- 
dius r* is large (r* ~ 10 6 cm. for neutron stars) which 
means that mor* = r*/A c ^> 1. 

Let us estimate the number of mor* for a neutron star 
using the following relation, 



Pc Pc 



(21) 



where P c is the pressure at the center of the neutron star 
and the equation of state P c = ujp c is assumed. 

Using the approximation Rq ~ p c and the constant 
equation of state parameter, u> ~ 0.1, we can rewrite r* 

— 1/2 

in the term of Hubble length r H = c/ Hq ~ R c as 



1/2 



,-i/a (RcY 



r*~i*n^=(^l r "~~(^) 12 ^o. ( 22 ) 



The faster shrinking of A c via Chameleon effect can be 
seen using Eq. (fT5|) 



A c = l/m 



Pa 



n+1 



(23) 



which means that the size of mor* depends on the density 
contrast in the following way 

/ v ra+1/2 

m r, = r*/A c ~ \^-\ ~ 10 43 ("+V2). (24 ) 



With the minimum requirement n > 0.9 for local gravity 
constraints, the growing mode at r = r* becomes 



- exp[10' 



(ill! 



(25) 



If C\ ^ or the growing mode is allowed, the initial 
condition (5?/>i(0) must be fine tuned to a fantastic level 
in order to compensate the enormous factor exp[10 60 ]! 

As we have shown, this tuning problem arises from the 
same criterion as Dolgov-Kawasaki time instability. The 
allowance of the incredibly small set of the initial condi- 
tion corresponding to the correct boundary condition of 
GR-likc solution should be considered as a serious theo- 
retical problem. The system which is highly sensitive to 
the initial value to the level of O(10 60 exp[10~ 60 ]) should 
not be considered as satisfactory. 

Let us note that setting C\ = from the beginning is 
not permissible because the continuity of the gradient of 
the solution ij){r) = ipi(r) + Sipi(r) at the center of the 
star needs 

V4(0) + Si'[(0) = ^(0) + tf^i(0) = 0, (26) 

which leads to 6ip[(r) = or C\ = —Ci whereas the 
for GR-like solution to hold, Vi(0) = V'o(0) = and 
mo ~ const are assumed. Then the perturbation 5ipi(r) 
can be rewritten as 



sinh(mor) 



mor 



(27) 



which gives the evolution of perturbation from centre to 
the surface of the star. We can estimate how the initial 
perturbation,5?/>i(0), at the center enhances as we move 
to the surface of the star, 



<fy>i(r*) « £Vi(0) 



601 



cxp[10 
2 x 10 60 



(28) 



The situation becomes worse for the growth of perturba- 
tions in case n > 2 with mor* ~ 10 107 . 

There are two ways to deviate from the first order it- 
eration solution tpi(0) depending on the initial sign of 
6tpi(0). For the positive <5?/>i(0), the perturbation drives 
i/>(r) toward curvature singularity ip{r) — 1 while in the 
case of (0) < 0, one can easily get ip(r) which is in- 
consistent with the observational constraints. Since ip{ r ) 
and i?(r) are in one to one correspondence by definition, 
the GR-like solution ip(r) cannot much depart from the 
value ipo (r) and tpi (r) for the entire stellar radius other- 
wise many known constraints of general relativity such 
as the one coming from double pulsar tests would not be 
satisfied. 

The simplest way to satisfy stringent test of general rel- 
ativity is provided by taking the Schwarzschild dc-Sittcr 
solution with small cosmological constant, 2XR C . From 
the 1 st iteration, the small interval of allowed deviation 
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V>i(r) - ip (r) 



/ D \ 2n+2 

4nA 2 (2n+l)(^J , (29) 



2n+2 



(30) 



where we have assumed n, A ~ 0(1). 

In case dipi(0) < 0, the maximum allowed deviation at 
r* can be approximated by 

2n+2 



<tyi(0) 



cxp [moT 
2m r, 



> _ 



(31) 



This requirement for GR-like solutions hold leads to an 
extreme fine tuning, for example, in case n = 0.9 (mo?'* = 
10 60 ), we need |^i(0)| < 6.3 x 10~ 104 exp[-10 60 ]. With 
this huge size of the scalaron mass inside neutron stars, 
the GR-like solution is highly unstable. Any small per- 
turbation from the exact solution, no matter how small in 
the physical sense, can cause a catastrophic effect such as 
the divergence of R or the solution cannot be compatible 
with observations Hit . 



Around the Big Ban 

n ~, m30 

I BBN 



nucleosynthesis (BBN), the density 
IO^'pa |22| giving rise to m^, ~ 10~ 4 Mp for 
n = 0.9 which is heavier than the typical inflaton mass 
10~ 6 Mp. The scalaron mass exceeds the Planck mass by 

29 



the factor 10 in case n = 2. 

In stars with nuclear matter density, p c ~ 10 p A , rn^p 
becomes 10 20 - 7 M P for n = 0.9 and 10 68 M P for n = 2. 
Hence, the length scales A c corresponding to this un- 
physical heavy masses are definitely shorter than Planck 
length , cquivalently, the gigantic numbers, m r* = r»/A c 
is also un-physical due to the uncontrollable Chameleon 
effect. 

The simplest way to remedy the fine tuning problem is 
provided by putting a cut off on by carefully chosen 
the maximum value of my, such mof* ~ 0(1) which can 
be achieved by adding (u/R c )R 2 term into the model 
under consideration (El 1 1 0| . 

Indeed, in the limit (R 3> R c ) in context with the 
Starobinsky model, we can use the approximation 



1 



A,_R_R &,RR 



- R 



(35) 
(36) 



IV. SCALARON MASS CUT OFF INDUCED 
BY (p/R c )R 2 TERM 

The fine tuning problem in f(R) theories is closely as- 
sociated with the large mass, the scalaron acquires in the 
high curvature regime which is essential for local grav- 
ity constraints to be evaded. It is, however, unsatisfac- 
tory that the scalaron mass can easily exceed the Planck 
mass[l|. Using the dimensional arguments, we can esti- 
mate the stellar radius of neutron star from the nuclear 
matter density, cosmological density and the Hubble ra- 
dius, 



r* 

y 



1/2 



r* = 10~ 43/2 x 1.3 x 10 26 m. (32) 



which gives the correct order of magnitude for the radius 
of neutron star, 0(1O 4 ) m. It is clear from Eq. (TTS)) that 
mass scale (length scale) corresponding to the scalaron 
mass increases (decreases) faster with density for any n > 
0, 



my, 




(33) 



As pointed out in Ref. ||, the scalaron mass can exceed 
the Planck mass even in the regime of the density of 
classical relativity. Let us estimate \fR c in the unit of 

M„, 



where A rr = dA ^/dR. The additional term (/// R C )R 2 
can provide a cut off on ~ {R c /(6n)y/ 2 when \x is 
chosen to satisfy the condition Q 



> fJ, > 



2(n+l) 



(37) 



In case of neutron stars, the upper limit /i <C p A /p c ~ 
10~ 43 is equivalent to (p/R c )R 2 <C R or the correction 
term should be small compared to the background curva- 
ture R. It should be noted that the lower limit is always 
satisfied for generic values of n. 

For the numerical calculation, we need not to stick to 
approximation (|3"T|) as can be calculated directly. As 
demonstrated in Ref.JT^], the carefully chosen p corre- 
sponding to the cut off on mass about the scale of neu- 
tron star density, can give rise to GR-like solution in 
for neutron stars. We may explain this observation by 
consider the heuristic argument (|32|) which gives r* ~ 

10 -43 / 2 i? c and the corresponding Compton wave- 
length for the mass cut off, A c ~ 1/my, ~ (6/i) 1 / 2 i? c 
such that 



— = \c 

m 



when 6p ~ 10 



-43 



(38) 



which clearly allows to avoid the catastrophic fine tuning 
problem in compact objects like neutron stars. 



VRc 



H 2.13/i x lQ~ i2 GeV/h 



1.22 x 10 19 GeV/c 



10 



-61 



(34) 



From Eq.Q and Eq. ([gg]). we find that > M p 
when p > 10 32 p A for n — 0.9 and p > 10 20 p A for n = 2. 



V. CONCLUSIONS 

We have examined the scalaron dynamics in the frame 
work of Starobinsky model. The curvature singularity is 
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generic to this class of models if they are to be consis- 
tent with local gravity constraints. The finite potential 
barrier between the de-Sitter minimum and the curva- 
ture singularity is a serious threat to models with dis- 
appearing cosmological constant. The problem becomes 
grave in high curvature regime in compact objects like 
neutron stars. In this case, the scalaron mass becomes 
larger than the Planck mass due to chameleon mecha- 
nism necessary for local gravity constraints to be evaded. 
This in turn becomes the root cause of instability prob- 
lem in the scenario under consideration. While evolving 
the scalaron from the centre to the surface along the ra- 
dius of neutron star, we need to stay close GR. Little 
perturbation of initial conditions at the centre can easily 
destroy the desired evolution. This condition for having 
GR-likc solutions requires extreme fine tuning of initial 
conditions, for instance, in case of n = 0.9, m^r* = 10 60 , 
we need \Sipi(0)\ < 6.3 x KT 107 exp[-10 60 ]. Situation 



further worsens for larger values of n. This intractable 
level of fine tuning throws a serious challenge to f(R) 
theories consistent with local gravity tests. 

The problem can be alleviated by introducing 
quadratic curvature term [iR 2 /Ro in the Starobinsky 
model which is equivalent to putting a cut off on scalaron 
mass corresponding to /i <C p A /p c ~ 10 -43 for nuclear 
matter density. This prescription, however, runs into 
problem if one asks for its compatibility with early uni- 
verse physics a la inflation which leads to much larger 
value of the scalaron mass. 
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